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Abstract 

This paper presents a new technique for location of transmission line faults using voltage and current 
measurements from one end of the faulted line. The method differs from past techniques in many respects. 
The main feature of the new technique is that it considers the influence of the remote-end infeed of the 
transmission line, the effect of the transmission line capacitance, fault resistance, prefault loading 
conditions, and the effect of mutual coupling between different phases of the line, and then computes with 
a microprocessor the distance to the fault point without any approximations. It uses recorded phase 
voltages and currents at the near end, then the fundamental components of the measured signals are 
extracted using microprocessor filters. The modal values of the extracted signals are computed then 
processed to indicate the precise location of the fault. 



Introduction 

Rapid power system repair and service restora- 
tion subsequent to a permanent fault on a trans- 
mission line are essential conditions of utility 
operation. With the advent of microprocessor 
based recording devices, utilities have shown in- 
creased interest in implementing accurate fault 
location techniques. Consequently, several meth- 
ods have been introduced and implemented dur- 
ing recent years. 

Takagi et al. [1], Eriksson et al. [2], Cook [3] 
and Wiszniewski [4] proposed techniques which 
use fault current distribution factors, and pre- 
fault and postfault currents and voltages from 
one line terminal; impedances of equivalent 
sources connected to the line terminals are used 
in the computation. In practice, from time to 
time the distribution factors are changed. Sev- 
eral approximations have also been introduced 
during the computation process. Richards and 
Tan [5] proposed an iterative technique which 
uses several system parameters for estimating 
the location of faults. They developed a sequence 
network approach that treats both the fault loca- 
tion and the fault resistance as unknown quanti- 
ties. Shunt reactance is not modeled, limiting the 
application of the algorithm to short lines. 
Schweitzer [6, 7] performed an iterative calcula- 



tion of the difference in phase angle between the 
total fault current and the remote-end feed cur- 
rent contribution. A sequence model of the line 
and remote-end source is used, but the mutual 
coupling between sequence component networks 
is ignored, thereby assuming that the phase 
impedances are balanced. Jeyasurya and Rah- 
man [8, 9] presented an algorithm which is based 
on using current and voltage information at both 
line ends and by solving the sequence network of 
the equivalent power system the location of the 
fault can be estimated. Johns et al. [ 10] proposed 
a technique which is based on modal analysis 
using voltage and current data available at both 
ends of the protected line section. Sachdev and 
Agarwal [11] used local digital impedance and 
relay current data as well as corresponding data 
from the remote end to calculate the distance to 
the fault point. This method gives considerable 
errors for certain locations and the fault resis- 
tance estimates are inaccurate for locations near 
the midpoint of the line where the fault current 
contributions from both line terminals are equal. 
Lawrence and Waser [12] developed an al- 
gorithm based on Z-transform techniques using 
time domain representation of voltage and cur- 
rent signals and the phase network model. 

This paper presents an accurate algorithm to 
estimate the location of a transmission line fault. 
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The method uses the fundamental frequency 
voltages and currents measured at one line ter- 
minal and takes into consideration the following 
factors: 

(1) remote-end current infeed; 

(2) exact transmission line representation 
(distributed parameter representation) is consid- 
ered in the computation to account for capac- 
itive currents; 

(3) prefault and postfault line currents and 
voltages are measured and stored and used for 
computation; 

(4) the contribution of current from relay loca- 
tion side to the total fault current at the fault 
location, assumed to be equal to an impedance 
ratio [2], is corrected; 

(5) no major approximations are carrried out 
during computation; 

(6) a modal analysis is used, permitting the 
inclusion of unbalanced system impedances 
caused by untransposed lines. 



Development of the algorithm 

The method to compute the distance to a fault 
point is based on the following equations in modal 
form [13,14]: 

Vp = Vs cosh(GZ) - Zg/s sinh(G/) 

/^p = {V'slZs)sinh{Gl) - /g cosh(GO 

where / is the distance to the fault point from the 
sending end of the line; Vg is the modal voltage at 
the sending end (postfault value), /g the modal 
current at the sending end (postfault value), Vp 
the modal voltage at the fault point (postfault 
value), /p the modal fault current (postfault 
value), V's the superimposed modal voltage at the 
sending end (modal voltage difference between 
prefault and postfault), and /g the superimposed 
modal current at the sending end (modal current 
difference between prefault and postfault), and 
7sF the modal fault current feeding the fault from 
the sending end side; Zg is the modal surge im- 
pedance, and Z the transmission line impedance 
per unit length; G = T + jP is the propagation 
constant. 

Referring to Fig. 1, for the postfault condition, 
the modal voltages and currents can be ex- 
pressed as 

Vp = Vs cosHGl) - /gZg sinh(G/) 

and the superimposed components of the modal 
voltages and currents can be expressed as 



"^S-end e = T*jP 
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Fig. 1. Single-line diagram of the power system: (a) prefault 
condition; (b) postfault condition (Gp = fault conductance): 
(c) superimposed component of fault current. 



/^p = (yg/Zg)sinh(G/) - /^ cosh(GO 
In matrix form we can write 



Vp = [cosh(G/)] Vg - [sinh(G/)]/g 


(1) 


Vgp= [sinh(G/)]YgV;; - [cosh(GO]/s 


(2) 


where 




Vp'=[Vp,,yp2, Vp,3] 




[cosh(G/)] =diag[cosh(Gi/),cosh(G2 0, 


cosh(Gg/)] 


Vg' = [Vs,,Vs2, Vg,,] 




[sinh(G/)] =diag[sinh(Gi/), sin^Gg/), 


sinMG;,/)] 


Zg = diag[Zsi , Zs2, Zgg] 




Yg = diag[yg,,yg„ Yg,,] 




••SF = I-'SFI > -'SF2> -' SFsJ 




^ S ~ 1-' SI I ■' S2 I 7 S3 J 





The voltage at the fault location, Vp, is related 
to the fault current /p by the fault constraint 
matrix Gp (where the superscript p indicates 
phasor values) as follows: 



^F 



GlVl 



The fault constraint matrices for different 
types of faults are indicated in Appendix A, 
whereas the distribution of the sequence compo- 
nents of currents at the fault location is indi- 
cated in Appendix B. 

Substituting eqn. (B-1) into eqns. (1) and (2), 
we get 

[sinh(G/)]YgVs- [cosh(G/)]/s 

= /eGp{[cosh(G/)]Vg-[sinh(GZ)]Zs/g} 
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[cosh(GO] {[tanh(G/)] Ys V's - VI's} 

= kGAcosHGl)] {VVs - [tanh(G/)]Zs/s} 

/2{Uys-[tanh(GZ)]Zs/s} 

= [GF-i]'{[tanh(G/)]Ysy^-U/^} 
Selecting one of the three modes (e.g. mode 1), 



( START ^ 



and letting 








F = [Gp.-i]' 




= 


/ll /12 /13 

/2I /22 /23 

_/31 732 /33. 






then 




Mysi--^sitanh(Gi07si] 
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L/ll 112/131 


"ta 

t8 
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inhCGiOYsiVsi- 
inh(G2/)ys2Vs2- 
inhCGa/) 733^83- 


T" ' 
T" 

-'82 

T" 



k[Vsr-Zs^tanhiGJ)hi] 

= AJtanh(GiZ)ysiVsi-/si] 
+ A2[tanh(G2/)ys2V^2-/s2] 
+ A3[tanh(G3/)ys2Vs3-/s3] 



(3) 



Equation (3) is a complex equation. It can be 
split into two separate equations which can be 
solved to get the two real unknowns k and I. 
Eliminating k, we get one equation in I, the 
distance to the fault location from the sending 
end, in the form 

/[tanh(Ti/), tan(Pi/), tanhiT^l), 

tan(PJ), tanhCTg/), tanCPgZ)] = 

This equation can be solved by the Newton- 
Raphson method for successive substitutions, 
where an initial approximation is assumed for /, 
that is, Iq, and an iterative process is carried out 
until the required convergence is achieved. This 
iterative process is governed by 

h^i = ik-f(h)inik) 

The main problem now lies in computing the 
matrix F. 

The flowchart of the fault location algorithm is 
given in Fig. 2. 



Test results 

The sample power system selected for the stud- 
ies is a 500 kV, 236 km line section. The basic 
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Fig. 2. Flowchart of the fault location algorithm. 

TABLE 1. Basic transmission line data 

No. of circuits 2 

No. of conductors per phase 3 

No. of earth wires 2 

Conductor position symmetry None 

Conductor resistivity (f) m) 3.2 x 10~* 

Earth wire resistivity (D) 2.69 x 10" 

Conductor strand diameter (cm) 0.32 

Earth wire strand diameter (cm) 0.32 

Outer diameter of earth wire (cm) 2.86 

No. of effective strands in phase conductor 54 

No. of effective strands in earth wire 54 

Earth resistivity (D m) 20 



transmission line data are given in Table 1. Two 
equivalent power systems were connected to 
buses R and S. Three types of fault were consid- 
ered (phase to ground, phase to phase, and dou- 
ble phase to ground). Faults were applied at 
different locations on the line under different 
operating conditions. The voltage and current 
phasors were calculated by the program, then the 
equivalent modal values were computed. Mode 2 
was used in the study. These data were used by 
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TABLE 2. Results of the study 



Fault type 


Assumed fault 


Calculated fault 


System 




location 


(p.u.) 


location (p.u.) 


condition 


Phase to earth 


0.6 




0.6 


Unloaded 


Phase to phase 


0.6 




0.6 


R|., = 50 n 


Double phase to earth 


0.6 




0.6006 




Phase to phase 


0.7 




0.7 


Light loading 


Phase to phase 


0.7 




0.7 


flp = ioa 


Double phase to earth 


0.7 




0.7005 


Power flow S to R 


Phase to earth 


0.9 




0.9 


Light loading 


Phase to phase 


0.9 




0.9 


i?H. = 20 n 


Double phase to earth 


0.9 




0.8999 


Power flow R to S 


Phase to earth 


0.2 




0.2 


Heavy loading 


Phase to phase 


0.2 




0.2 


Rf = Q 


Double phase to earth 


0.2 




0.1998 


I'ower flow S to R 



the fault locating algorithm to evaluate the dis- 
tance to the fault from the sending end. Table 2 
summarizes the results of this study. 



Conclusion 

This paper has presented an accurate al- 
gorithm for locating transmission line faults. The 
algorithm has been tested for a variety of simu- 
lated fault conditions, including double-circuit 
mutual coupling effects, various remote source 
equivalent conditions, various prefault steady- 
state conditions, and in the presence of apprecia- 
ble fault resistance. These simulation studies 
indicate that the proposed method is accurate for 
all types of shunt fault at different locations 
under different operating conditions. The maxi- 
mum error was found to be 0.1%. 



Appendix A 

Fault constraint matrices 

Referring to Fig. A-1, the fault constraint 
matrix Gp is given by 

(a) for a line to ground fault 





"1 0" 






Gl = y, 



.0 0_ 




(b) for a line to line fault 




ro 01 




dP _yF 



-1 


-1 
1 





D^g 



Rr 



L-G fault 




2-L-G foult 



R^ 



L-L fdull 



pRrg^FQ ^F- 



3-L fault 



tj^F n^F 



3-L-G fault 



Fig. A-1. Types of fault. R^ is the arc resistance from phase to 
tower plus the effective tower footing resistance for a single line 
to ground fault, and the effective tower footing resistance for 
other types of ground fault. R^ is the arc resistance between 
faulted phases. 



(c) for a three-phase fault 



^^ 3 



2 -1 -1 
-1 2 -1 
-1 -1 2 



(d) for a double line to ground fault 
1 



Gl = 



Rp^ + 2RfRg 















R,+Rg 


-a. 





-n. 


R, + Rg_ 



(e) for a three-phase to ground fault 

yo + 2yp Jo-Jf Jo~Jf n 

Gi = - jo-Jf yo + ^yp yo-yp 

yo - Jf Jo - Jk 3*0 + 2Jf 
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where 

y, = lliRr + 3R,) 

The phase quantities F" and /p are related to 
the modal quantities V and / by the eigenvector 

matrices Ci and C2: 

P = cj 

Since we have 

where S^ is a real quantity representing the fault 
resistance, 

Iyig for a line to ground fault 
jp/2 for a line to line fault 
jp/3 for a three-phase fault 

For a double line to ground fault, 

S, = l/(i?p2 + 2R^^RJ 



Gp = Rf 



In this paper Rg is taken to be equal to Rp. 
For a three-phase to ground fault, 



' SFa^" SFao*'sFa1 *'sFa2 



0" 




1 


+ Rg 


1 
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-1 





-1 


1 



'Si=Jf 



Gl = 



Yq-Jf 

Syr 



1 1 1" 
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1 1 1 





1 0' 
10 
1 



Appendix B 

Distribution of sequence components of currents at 
fault location 

Referring to Fig. B-1, the phase relation be- 
tween the fault current /p and the fault current 
from the sending end is governed by the current 
distribution factors k^, k^, and /jg of the sequence 
components: 

^0 — •'sfoZ-'fo 
f^i — •'sfiZ-'fi 

1^2 ~ ■'SF2/-'F2 

We can assume to a high degree of accuracy that 
kj = kz, so that 




'sFb='sFbo*'sFb1*'sFb2 



4J'Fb='Fo^lFl- 



'F2 



L-L fault 



'sFa = 's&ol'sFat *'sFa2 



'5Fb='SFbo*'SFb1*'SFb2 



MK 



Fao* 
'Fal* 
lFa2 



3-L fault 




'sFa='SFao*'SFa1*ISFa2 



lFa2 



3-L-G fault 

Fig. B-1. Distribution of sequence components of currents at the 
fault location. 



■'SfZ-'F ~ (-'SFO + -'SFl + -'sF2)/(-'F0 + ■'FI + ■'F2) 
= L^o-'fO + f^l\^Fl + •'F2)J /■'f 
= «1 + {Rq — /ji)ipo/iF 

(a) For a single line to ground fault, 

ip = oi Fo 

and 

(b) For a line to line fault, 
/fo = 

■'sFb/-'Fb = Ki = «LL 

(c) For a three-phase fault, 
/fo = 

and 

-'SFaMFa = «1 = «3L 

(d) For double line to ground and three-phase 
to ground faults, 

■'SFaMFa = «1 + («0 ~ ^1 KpOa /-'Fa 
-'SFb/-'Fb = «i + («o ~~ «l)-'FOb/-^Fb 
I'sFclhc = ^1 + (^0 - ki)Iyoclhc 

But 

^FOa — ^FOb — -'fOc ~ -'fO 
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and 



T" 



■ ■'SFOb — ''sFOc ~ ^O-'fO 



Let 

{ko-k,)rspoiko = Co 

then 



T" 

J- SFa " 


-Co" 




'^Fa" 


T" 

J- SFb " 


-Co 


= k. 


Iv, 


T" - 


-Co 




Jv._ 



or 



IsFsL 




'W 




"1" 


Isvh 


= k. 


h. 


+ Co 


1 


T" 




hr_ 




1 



In the case of double line to ground and three- 
phase line to ground faults, k^ is assumed to be 
equal to kg, that is, Cg = 0. 

Hence, /g^ can be represented by a real frac- 
tion, ky, of the total fault current /p (e.g. ki^a, 
^LL. ^3L> ^2LGi ^3lg)> that Is, the angle of the fault 
current /p and the fault current from the sending 
end, /gp, are assumed to be equal. 

In phasor form 



••SF 



^ k,Jl 



■* SF — '~'l '^f'-'f ' F 

In modal form 

/gp = 0|«pC2 (jpCi Vp 

= A:GpVp 
where 

Gp = Ca^'GpC, 

= c,'Gs;x, 



(B-l) 
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